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Abstract 

A new method for the preliminary design of con- 
trolled space structures is presented. The method 
coordinates standard finite-element structural anal- 
ysis, multivariable controls, and nonlinear program- 
ming codes and allows simultaneous optimization of 
the structures and control systems of a spacecraft. 
Global sensitivity equations are a key feature of this 
method. 

The preliminary design of a generic geostationary 
platform is used to demonstrate the multidisciplinary 
optimization method. Fifteen design variables are 
used to optimize truss-member sizes and feedback- 
gain values. The goal is to reduce the total mass of 
the structure and the vibration control system while 
satisfying constraints on vibration decay rate. Incor- 
porating the nonnegligible mass of actuators causes 
an essential coupling between structural design vari- 
ables and control design variables. 

The solution of the demonstration problem is 
an important step toward a comprehensive controls- 
structures integrated design capability. Use of global 
sensitivity equations helps solve optimization prob- 
lems that have a large number of design variables 
and a high degree of coupling between disciplines. 

Introduction 

Future NASA missions will include large space 
structures with control systems to damp out vibra- 
tions excited by pointing maneuvers. Preliminary de- 
sign of these spacecraft is complicated by a high de- 
gree of coupling between the control and structural 
analyses. Specifically, changes in the structure im- 
pact the control system design by modifying both 
the plant to be controlled and the expected exci- 
tation. At the same time, changes in the control 
system impact the structural design by modifying 
the number, mass, and location of actuators. The 
preliminary design of controlled space structures is 
one aspect of a much broader program called the 
NASA Controls-Structures Interaction (CSI) Tech- 
nology Program (ref. 1). 

The CSI program encompasses a variety of 
research areas. The primary objectives of the pro- 
gram are (1) to develop analytical methods that char- 
acterize the performance of controlled space struc- 
tures, (2) to develop controls-structures integrated 
design (CSID) methods, and (3) to improve ground 
test methods to better predict on-orbit system per- 
formance and validate the analytical methods. The 
first objective not only emphasizes the interaction 
between the structure and various control systems 


but also quantifies the performance benefits and op- 
erational restrictions that result from this interac- 
tion (refs. 2 and 3). The second objective exploits 
advances in structural analysis, multivariable con- 
trol, and multidisciplinary optimization for the pre- 
liminary design of large, flexible spacecraft (refs. 4 
to 8). The third objective emphasizes hardware im- 
plementation and determines the validity of analyt- 
ical assumptions. Thus, this objective requires new 
methods for accurate static and dynamic modeling 
of structural and control system hardware, innova- 
tive and physically realizable control strategies, and 
new methods for experimental system identification 
and verification (refs. 2, 3, and 9). 

Research in CSID has grown during the past 
10 years. Initial research studied the mechanisms 
that link structural characteristics to controlled per- 
formance by using simplified models and control the- 
ory. Reference 4 contains a survey of this work 
and concludes that strongly coupled CSID problems 
need attention. Recent research considered prelim- 
inary design problems that are complex enough to 
capture the most important characteristics of de- 
sign problems associated with actual mission hard- 
ware. For example, researchers at NASA Langley 
Research Center have developed mathematical opti- 
mization procedures that successfully reduce surface 
distortion errors for large space antennas (ref. 5), tai- 
lor structures and control systems for reduced power 
consumption (ref. 6), and improve the fine-pointing 
performance of large space platforms while reducing 
mass (refs. 7 and 8). 

The research reported in this paper addresses 
problems for which there is implicit coupling between 
structural design variables and control design vari- 
ables. These optimization problems are challenging 
because they involve eigenvalue solvers and transient 
response analyses, which can be computationally ex- 
pensive. Moreover, these analyses must be iterated 
until all structural and control response quantities 
converge. 

The current research is unique because it coordi- 
nates commercial structural analysis codes, 
multivariable control codes, and optimization codes 
using the UNIX command language. This method 
can also be adapted for a variety of optimization 
problems associated with controlled space struc- 
tures. The optimization method builds on recent ad- 
vances in multidisciplinary analysis and optimization 
(refs. 10 to 12). Two key features of this research are 
the division of the system engineering problem into 
subproblems and the use of general-purpose anal- 
ysis and optimization computer software for pre- 
liminary design studies. Because the disciplinary 



analyses are tightly coupled, the global sensitivity 
equations (GSE) are used to calculate global deriva- 
tives of the response quantities with respect to the 
design variables. (See ref. 10.) For completeness, the 
GSE approach is compared with a more conventional 
approach. 

The example used to illustrate the optimization 
method is the preliminary design of the structure 
and a vibration suppression controller for a geosta- 
tionary platform subjected to slewing maneuvers. In 
this problem, truss-member sizes and controller gains 
are chosen to minimize spacecraft mass, which in- 
cludes vibration control actuator mass. The mass 
of the vibration control actuators is a parameter 
that, influences the spacecraft dynamic characteris- 
tics and a control system parameter that determines 
the maximum torque available for vibration suppres- 
sion. Thus, the structural and control analyses are 
coupled by the presence of the nonnegligible actuator 
mass. 

This paper develops and demonstrates a new 
technique for CSID. The design procedure is sum- 
marized, and alternative approaches for calculating 
global derivatives are discussed. The reference con- 
figuration and the optimization problem are pre- 
sented with a brief discussion of the coupled disci- 
plinary analyses. Typical results are presented, and 
a comparison between the GSE approach and a con- 
ventional approach is given. 

Symbols and Abbreviations 

CEC collocated elastic control 

CPU central processing unit 

CSI controls-structures interaction 

CSID controls-structures integrated design 

EAL Engineering Analysis Language 

G p , G r position- and rate-gain matrices 

used to define collocated elastic 
control law 

GSE global sensitivity equations 

g vector of 12 design variables used to 

define Gp and G r matrices 

I/O input and output 

J spacecraft inertia matrix, kg-m 2 

Afact total mass of vibration suppression 

actuators, kg 

A/ s total mass of truss structure, kg 


n number of modes used in reduced- 

order model of spacecraft 

obj objective function 

r vector of three design variables used 

to define radii of truss-element cross 
section, m 

v vector of design variable values 

associated with a specific stage in 
optimization procedure 

x arbitrary design variable, x £ {r.g} 

y arbitrary constraint function 

a distance between points in design 

space 

<5 required vibration decay rate 

A 2 n x 1 vector of closed- loop eigen- 

values, rad 2 /sec 2 

& 3 x n mode-slope matrix that 

contains rotational components 
of structural eigenvectors; super- 
script 0 indicates maneuver actu- 
ator location and superscripts 1 
and 2 indicate first and second 
vibration-suppression-actuator 
locations 

cj 2 n x 1 vector of natural eigenvalues 

of the structural model, rad 2 /sec 2 

Multidisciplinary Optimization 
Procedure 

CSID Method 

There are three phases in the CSID process: anal- 
ys\s, calculation of global derivatives, and optimiza- 
tion. Analysis involves iteration between structural 
and control analyses until output quantities, such as 
structural mass and closed-loop response, converge. 
The global sensitivity calculation produces deriva- 
tives of response quantities for each design variable. 
The optimization phase searches for new values of 
design variables that reduce the objective function, 
reduce the degree of constraint violation, or both. 
Since the coupled analysis is expensive, linear ap- 
proximations to the objective and constraint func- 
tions arc used in the optimization phase. Limits on 
the allowable change in design variables restrict the 
search algorithm to a domain in which the linear ap- 
proximations are appropriate. When no furl her re- 
duction of the approximate objective function can 
be achieved without constraint violation, the coupled 
analysis and the global derivative calculations are re- 
peated for the next set of design variables, and the 
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process continues. One repetition of analysis, deriva- 
tive calculation, and optimization is referred to as 
one cycle. The CSID process continues until a pre- 
scribed maximum number of cycles is readied or until 
all constraints are satisfied and the change in the ob- 
jective function is less than a prescribed convergence 
tolerance for two successive cycles. 

A conventional three-phase CSID process is il- 
lustrated schematically in figure 1(a). The coupled 
analysis is represented by an iteration between the 
controls and structures computer codes. This cou- 
pled analysis is performed once at the nominal design 
point and is repeated after each design variable is 
perturbed. The global derivatives are approximated 
with a forward-difference formula. The optimization 
phase is represented by an iteration between the lin- 
ear approximation and optimization computer codes. 
The linear approximation computer code calculates 
approximate values of both the objective function 
and the constraints by a first-order Taylor series ex- 
trapolation from the nominal solution by using gra- 
dient information provided by the global derivative 
phase. The optimization computer code is a general- 
purpose optimization program. 

Figure 1(b) illustrates an alternative CSID ap- 
pioach in which the coupled analysis is used once per 
cycle. The controls and structures computer codes 
are modified to calculate the partial derivatives of 
output quantities with respect to input quantities. 
The 1 global derivatives of response quantities with re- 
spect to each design variable are calculated with this 
local sensitivity derivative information. 

1 he calculation of global derivatives in the alter- 
native approach discussed in appendix A is imple- 
mented through the use of global sensitivity equa- 
tions (GSE) described in references 10 and 13. The 
global derivatives are expressed as the solution to 
a set of linear equations whose coefficients are local 
sens i t i vi ty derivat i ves . 

To calculate all global derivatives, it is necessary 
to solve a series of GSE problems that have the form 
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where I is the identity matrix, S is the set of all 
structural outputs, C is the set of all control out- 
puts, and x is any global design variable. Standard 
mathematical computer software can be used to solve 
equation (1) for each design variable x. The multi- 
ple solutions arc calculated efficiently since only the 
right-hand side changes. 


In this paper, both the conventional and the GSE 
approach are used. The conventional approach is rec- 
ommended for uncoupled problems. In this ease, the 
disciplinary analyses can be computed without iter- 
ation, and finite-difference approximations to global 
derivatives are tractable. The conventional approach 
is also recommended for coupled problems for which 
the number of structural modes used in the control 
analysis is large' compared with the number of design 
variables. In this cast', the dimensions of the GSE! 
matrix can become large, and the cost of evaluating 
the local derivatives can become excessive. On the 
other hand, the GSE approach is ideal when the num- 
ber of design variables is so large that repeated per- 
turbations of the coupled analysis is computationally 
prohibitive. The size of the GSE matrix is unaffected 
by the number of design variables x. Increasing the 
number of design variables means that equation (1) 
must, be solved for additional right-hand sides. This 
impacts the computational cost only when the local 
derivatives in the added column vectors are expensive 
to calculate. 

Demonstration Problem 

The preliminary design of a generic geostation- 
ary platform is used to develop and demonstrate 
the CSID methods. The goal is to reduce the total 
mass ol the structure and the vibration control sys- 
tem while satisfying constraints on vibration decay 
rate. Standard finite-element analysis, multivariable 
control, and mathematical programming routines are 
coordinated by a new multidisciplinary optimization 
scheme. 

The reference configuration, shown in figure 2, 
represents a large geostationary platform, which con- 
sists of two antennas connected by a bus structure. 
The platform has three sets of actuators; each set 
supplies three-axis torques. Two of the actuator 
sets work together for vibration suppression control. 
These sets are mounted at the intersection of the cen- 
terline of the bus structure and the axis of symmetry 
ol each antenna. The third set of actuators is used 
to perform rigid-body rotational maneuvers and is lo- 
cated at an arbitrary joint in the bus structure near 
the center of mass. A finite-element, model is assem- 
bled from beam elements to represent a graphite- 
epoxy truss structure. Three groups of beam ele- 
ments arc defined: the bus, the antennas, and the 
antenna supports. All beam elements have the same 
wall thickness, but the outer radius may vary from 
one group of elements to another. 

To obtain a tractable integrated design problem, 
certain simplifying assumptions are made. The loca- 
tions of the three actuator sets are fixed, the mass 
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of the maneuver actuator is ignored, and the inertia 
matrix J is calculated at the fixed maneuver loca- 
tion rather than at the variable eenter-of-mass lo- 
cation. Vibration-suppression-actuator masses are 
included in the model, but their rotational inertias 
are not included. The masses of scientific payloads, 
power-generating and power-distribution hardware, 
computers, and other supporting equipment are ig- 
nored. The mass and dynamic characteristic s of the 
joints between truss elements are also neglected. 1 he 
stiffness contribution of the surface material of the 
antenna reflector is neglected, but the mass is mod- 
eled approximately. Minimum gage constraints on 
the truss-member sizes arc 1 also used instead of stress 
and buckling constraints. The structural model is 
discussed in appendix B. 

To complete the structural model, the masses of 
the two vibration suppression actuators must, be de- 
termined and incorporated into the model. A con- 
straint in the control law design requires that the vi- 
bration suppression actuators to have the same mass. 
To determine the actuator mass, first a torque his- 
tory characteristic of a minimum-time rigid-body ro- 
tational maneuver is applied to the elastic spacecraft 
with the vibration suppression loop closed. The mass 
of each act viator is calculated based on the sum of the 
maximum torques about each axis required to sup- 
press vibrations that occur after this reference ma- 
neuver. The total mass is obtained by multiplying 
the number of actuators (two) by the mass per actu- 
ator. The actuator mass calculation is discussed in 
appendix C. 

Controls-structures integrated design of the ref- 
erence configuration is implemented with a mathe- 
matical programming method used to adjust the val- 
ues of 15 design variables. Three design variables r; 
control the outer cross-sectional radius of the gioups 
of beam elements that comprise the bus structure 
(i = 1), the antennas (i = 2), and the antenna sup- 
ports (i = 3). The remaining 12 design variables 
g x (i = 1,2, .... 12) uniquely determine the elements 
of the position- and rate-gain matrices for the vibia- 
tion control system. 

The optimization problem is to minimize the total 
mass of the platform while satisfying vibration decay 
requirements. This problem can be stated in terms 
of the total mass of the truss structure Af s , the total 
mass of the vibration suppression actuators A/ ac t , 
and the complex- valued, closed-loop eigenvalues X t . 
Minimize 

A/,s -b A /^ t (^ a ) 

subject to 

Rc(A i)<6 (i = l,2,...,2n) (2b) 


where t) is the required decay rate and n is the 
number of modes used in the reduced-order model of 
the spacecraft. The calculation of the real part of the 
closed-loop eigenvalues is discussed in appendix C. In 
the remainder of the paper, \ t is used to stand for 
Re(A / ). 

This design problem is particularly challeng- 
ing because the structural and control analyses are 
tightly coupled. Changes to the structural design 
variables (truss sizes) change the open-loop plant and 
thus the closed-loop response to the reference ma- 
neuver. On the other hand, changes to the control 
design variables (the gains) determine the vibration- 
suppression-actuator torques, which are linearly re- 
lated to the actuator masses used in the structural 
model. 

Controls-Structures Analyses 

This section briefly describes the structural and 
control analyses and explains the successive substi- 
tution algorithm used to solve the coupled problem. 
The inputs to the coupled controls-structures analy- 
ses are the design variables, and the outputs are A/ s , 
A/ art . and A used in equations (2). 

The structural analysis consists of finite-element 
analysis of the spacecraft structure through the use of 
the Engineering Analysis Language (EAL) computer 
code. (See ref. 14.) The mass and stiffness matrices 
are assembled based on the current values of the 
design variables v t and the mass of the actuators 
A/ ac t . A real symmetric eigenvalue problem is solved 
to obtain the characteristic modes and frequencies 
of the truss structure. Outputs from the structural 
analysis include the first n eigenvalues and the 
first n mode slopes (i.c., rotational components of the 
eigenvectors) at the maneuver actuator location ^ 
and at the vibration-suppression-actuator locations 
<j>0) and In addition, the mass of the bare 

structure A/ s and the inertia matrix J associated with 
the maneuver actuator location are output. All these 
outputs, except for A/s, are inputs to the control 
analysis. 

The control analysis determines the transient 
closed- loop response of the spacecraft to a reference 
maneuver and searches for the maximum vibration 
control torques and the times at which the\ occur. 
The maneuver is assumed to be linear since it rotates 
the spacecraft through small angles during a long pe- 
riod of time. Thus, the nonlinear coupling of rigid 
and elastic motion that is typical of rotational ma- 
neuvers is not modeled. The maneuver represents a 
typical disturbance that might be encountered by the 
geostationary platform. The vibration suppression 
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system employs collocated elastic control (CEC) with 
a single pair of actuators constrained so that their net 
torque output is zero. The robust dissipative control 
law used by CEC does not affect rigid-body motion 
and guarantees stability despite unmodeled dynam- 
ics and parameter uncertainty (refs. 15 to 17). The 
torque required by the CEC actuators is a function 
of the angular deformations and angular deforma- 
tion rates at the actuator locations and the values of 
the 12 design variables (j , , which uniquely determine 
the position-gain matrix and rate-gain matrix 
G, The complete control analysis is contained in 
appendix C. 


The actuators are sized based on the peak torques 
required to suppress the elastic motion t hat remains 
after the reference maneuver. The mass of actuators 
needed to provide the peak torques can be interpo- 
lated from an empirical table of mass versus torque 
for available actuators. For the present study, how- 
ever, a linear relationship between actuator mass and 
maximum torque output is assumed. Outputs from 
the control analysis include A/ ;ir j and the 2 // complex 
closed-loop eigenvalues A. The actuator mass M iWl 
is an input to the structural analysis. 

The solution to the coupled controls-structures 
analysis problem is found by iteration. Figure 3 
shows a flow chart of the iterative procedure that 
starts with an initial estimate of the actuator mass. 
This value of A/ art is incorporated into the struc- 
tural model. The required actuator mass is deter- 
mined from the t ransient response calculated in the 
control analysis. If the value' of A/ art used in the 
structural model and the value 1 calculated by the 
control analysis are within 1 percent of each other, 
then the procedure terminates. Otherwise, the re- 
quired actuator mass is used in the structural analy- 
sis (i.e.. the mass matrix is updated) and the iteration 
continues. 


The above procedure was modified to improve ef- 
ficiency. 1 he computational cost of the procedure 
is dominated by the cost of solving the structural 
eigenvalue problem. Fortunately, calculating the sen- 
sitivity of the solution of the eigenvalue problem to 
changes in lumped masses (e.g., du> 2 /()M HVi ) is rel- 
atively inexpensive. (Sec 1 refs. 18 and 19.) Then- 
fore, the change in structural outputs with a change 
bi A /act be estimated when the change is suffi- 

ciently small. Currently, the full structural analysis 
is repeated only when the change in actuator mass 
exceeds 10 percent of the initial estimate. Otherwise, 
the structural outputs are estimated from sensitivity 
information. 


The modification for efficiency required one addi- 
tional change in the procedure. If the change' in A/ art 
is less than 1 percent but is based on an estimate' of 
the structural outputs, then one additional oxeeu- 
t-ion of the structural and control analyst's is needed 
to confirm convergence. During the' additional exe- 
cution of the structural analysis, the' loc i al derivatives 
(e.g., Dlo- / dr 2 ) of structural outputs with respect to 
design variables are' calculated, along with the' local 
derivative's with respect to A/ a , . 

GSE and Optimization 

This section contains a summary of the global 
derivative' calculation and the optimization process 
required to complete one cycle of the new CSID 
procedure'. These calculatiems. which take place 
afte'r the controls-structures analysis has converged, 
assume that value's of the output epiaut.it ies J. oj. 
#«». *d>, A/ s , Mu, - and A are available. 

To solve the GSE problem (cep (])). all local 
derivative's that are' the 1 e*o('ffie‘i('iits of the matrix 
equations must be evaluated. That is. feu oae'h e*on- 
tributing analysis, the partial derivatives of each out- 
put quantity with respect, to e'aeh input quantity 
must be calculated or estimated. For the demon- 
stration problem, the contributing analyses are eon- 
trolis and structure's, and the input and output quan- 
titie's are 1 shown in figure 1 4. All partial derivatives 
a.i e seeded according to the 1 scheme suggested in rof- 
erence 20 and placed in the GSE as discussed in 
appendix A. 

I he sot of linear GSE (exp (1)) is assembled and 
solved to determine the' global derivatives of A/ s . 
AAirt «. and A wit h respect to e'aeh design variable 1 . 
For example, one of these 1 linear expiations is 


<IM u- t . 

dr*) 


n 


= £ 


()Af iir i duJ, 
Oujj dr 2 


y- ^A/ a ( t d J t 
OJi dr 2 


2 

EE 

/=i./=o o<i>j 


U) 


d.0 


U) dr 9 


( 3 ) 


Equation (3) reflects the fact that the required ac- 
tuator mass is not influenced directly by the truss 
sizing variable r 2 (i.e., 0M iici /dr 2 = ()) but is infhi- 
encxxl indirectly bv the eigenvalue's and the elements 
of the inertia and moeie-slope matrices. Tlie notation 
df dr 2 signifies the gle>ba.l derivative with respect to 
r 2* Thus, the global derivatives quantify the effect 
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of coupling between the local control and structural 
analyses. 

Once the global derivatives are available, they 
can be used to provide linear approximations to 
the objective function and the normalized constraint 
function y t . For example, given 

obj = M s + M a( -t 

y, = 1 - ^ (* = 1,2 2n) 

where y t < 0 is a feasible constraint, and the required 
decay rate' b is a negative real number; therefore, the 
effect of changing an arbitrary design variable x by 
the amount Ax is approximated by 



obj (a: + Ax) ~ M s 


dM H 

dx 


Ax 


, r dA/jict A 
-f A/jict H : Ax 


Vii x + Ax) « 1 


dx 

X t dX Ax 

b dx b 


( 5 ) 


The solution of the optimization problem in 
equations (2) is accomplished by linking the linear 
approximation routine with the standard nonlinear 
programming code CONMIN (ref. 21). A nonlinear 
programming code was selected because future CSID 
projects may require nonlinear objective and con- 
straint functions. 

Results and Discussion 

Demonstration Overview 

All results in this section were generated with the 
GSE approach. The impact of choosing the GSE ap- 
proach over the conventional approach is discussed 
in the section “Comparison of GSE and Conven- 
tional Approach. The implementation details are 
discussed in appendix D. 

The proper operation of the optimization pro- 
cedure is demonstrated through convergence to the 
same solution from two different initial design points. 
In each case, the total mass of the reference config- 
uration is minimized with respect to 15 design vari- 
ables (eqs. (2)). The decay rate requirement is b = 
—0.03. The reference maneuver rotates the space- 
craft 20° about all axes simultaneously in 10 sec. 
Limits are set so that a design variable may change 
no more than 10 percent during a single cycle. The 
number of modes and frequencies that describe the 
open-loop model was n = 20, and 0.5 percent modal 
damping was assumed. 


Case 1: slightly infeasible initial design. 

Demonstration case 1 has an initial design that does 
not fulfill the decay rate requirement. The initial val- 
ues of the truss sizing variables are r\ = 4, r<i — 4, 
and r 3 = 8 cm. These sizes yield an initial struc- 
tural mass of 1158 kg. The frequencies of the first 
five elastic modes are 0.668, 0.977, 1.35, 1.48, and 
1.48 Hz. The fourth and fifth vibrational frequen- 
cies are identical (i.e., they are repeated eigenvalues) 
because of structural symmetries in the spacecraft 
antennas. Both closely spaced and repeated frequen- 
cies are characteristic features of large space struc- 
tures and are considered in this demonstration prob- 
lem. The initial values of the position- and rate-gain 
matrices arc 


5100 

800 

700 

800 

5000 

700 

700 

700 

4900 


( 6 ) 


The closed-loop response that results from these 
gains indicates an actuator mass requirement of only 
10.69 kg. The first five closed-loop eigenvalue pairs 
are -0.021 ± 4.2j, -0.031 ± 6.1j, -0.049 ± 8.5j, 
—0.047 ± 9.3/, and -0.047±9.3j. Only the first pair 
violates the constraint that Re(A) < —0.03. 

Figure 5 contains optimization histories for demon- 
stration case 1. Figure 5(a) shows the history of 
the objective (i.e., total mass) with respect to the 
optimization cycle. The CSID procedure increases 
the total mass during the first five cycles and de- 
creases the total mass at each cycle thereafter. The 
initial increase in total mass is associated with a 
rapid decrease in the constraint violation, as indi- 
cated in figure 5(b). This figure shows that the max- 
imum eigenvalue constraint max[Re(A)] decreases 
smoothly to the constraint boundary in five cycles. 
The maximum constraint value oscillates around the 
constraint boundary beginning at cycle 15. Such be- 
havior is common when a linear approximation to 
nonlinear constraints is used. 

Figures 5(c) to 5(c) contain histories of selected 
design variable values. Only the diagonal elements 
of the gain matrices G p and G r arc plotted. The 
off-diagonal elements never change significantly from 
their initial values. A comparison of these histories 
with those in figure 5(a) shows that the initial in- 
crease in total mass is associated with increases in 
all three structural design variables and in rate gain. 
The subsequent decrease in mass is associated with a 
decrease in the radii of the bus truss. A decrease in 
rate gain also contributes, but this decrease lags the 
structural changes. This effect implies that changes 
in the structure not only decrease mass but also 
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reduce the amount of control effort required to meet 
decay rate constraints. 

Figure 5(f) is similar to figure 5(a) but uses a 
different scale for the actuator mass. This figure 
indicates that the actuator mass requirements are 
changed by ±25 percent from cycle to cycle, even 
when the change in design variables is much smaller. 
This behavior emphasizes that the amount of con- 
trol effort required is a function of diverse factors, 
such as dynamic and inertial properties of the struc- 
ture plus adjustments to the control system. More 
importantly, figure 5(f) indicates that an increase in 
actuator mass from about 10 to 50 kg allows a re- 
duction of about 300 kg in total mass. This favor- 
able trade-off of structure mass for control effort is 
the major benefit of an integrated controls-structures 
optimization procedure. 

Case 2: highly infeasible initial design . 

Demonstration case 2 has an initial design that has a 
smaller total mass and a greater constraint violation 
than the first case. The only difference between the 
cases is that initial values of the truss sizing variables 
are each set to 2 cm. This causes the first seven pairs 
of closed-loop eigenvalues to violate the decay rate 
constraint . Figure 6 contains a convergence history of 
mass with respect to optimization cycle. This second 
case converges to the same total mass as the first case, 
namely, about 850 kg with an actuator mass of about 
50 kg. 

One way to compare the optimization histories 
of the two cases is to plot the maximum eigenvalue 
constraint versus total mass. (See fig. 7.) In this 
figure, movement downward and to the left signi- 
fies improvement of the design from one cycle to 
the next. The two demonstration problems start at 
different initial conditions and follow different con- 
vergence paths. However, the final conditions are 
similar. This similarity shows that the optimization 
procedure performs well (consistently produces an 
improved design) and that the final solution is pos- 
sibly a global minimum. 

Although the optimization procedure converges 
to equally satisfactory solutions from different initial 
design points, the convergence path of demonstration 
case 1 is less direct than that of case 2. The path 
of case 2 makes consistent progress toward the final 
solution, while the path of case 1 initially adds an 
unnecessary amount of mass. Two of the structural 
design variables increase during the early optimiza- 
tion cycles but return to their initial values at the 
end of the optimization process. (See fig. 5(c).) This 
observation, and the previously noted oscillation of 
the maximum constraint value around the constraint 


boundary, suggests that the quality of global deriva- 
tives should be investigated. The derivatives may be 
accurate, but the move limits allow design variable 
values to change so much that the linear approxima- 
tion is a poor representation of the coupled analysis. 
Alternately, the global derivatives calculated by the 
GSE approach may be deficient. 

A simple test can be used to determine the quality 
of global derivatives that were calculated in demon- 
stration case 1. Two design points vj and v *2 that 
correspond to successive optimization cycles are se- 
lected. A coupled analysis is performed at each de- 
sign point and at evenly spaced intermediate points 
so that: 

Vnow = av 2 + (l-tt)vi (n = 0, 0.1, .... 1 .0) (7) 

where v now is a sequence of intermediate points be- 
tween the chosen design points vj and V 2 - Fig- 
ure 8(a) compares the actuator mass calculated by 
the coupled analysis with the mass predicted with 
nominal values and global derivatives at the point 
V]. In this case, and are the design points as- 
sociated with cycles 18 and 19 in figure 5. The differ- 
ence between actual and approximate A/ act is small 
at n = 0.1, but this difference is large at o = 1.0. 
This change indicates that the global derivatives are 
calculated correctly but that the move limits need 
to be small if the linear approximations for actuator 
mass are to remain valid. Fortunately, this discrep- 
ancy is less severe for the prediction of closed-loop 
eigenvalues. Figure 8(b) contains the corresponding 
comparison of the real part of the maximum eigen- 
value with the linear approximation to this value. 
The behavior is more nearly linear, and the error 
caused by the current move limits is acceptable. 

Assessment of GSE Approach 

This paper investigates the use of the GSE ap- 
proach to provide global derivatives for CSI opti- 
mization problems. As previously discussed, the 
GSE approach provides high-quality derivative in- 
formation. In this section, the strengths and weak- 
nesses of GSE are discussed from an implementation 
standpoint. The strengths involve numerical, effi- 
ciency, human judgment, and organizational issues; 
the weaknesses include implementation complexity 
and computer memory requirements. 

The numerical argument for the GSE approach is 
that it is relatively insensitive to the degree of con- 
vergence of the iterative (coupled) analysis. The it- 
erative procedure (fig. 3) terminates based on the 
percentage change in the mass of actuators. The 
current study experimented with 1 and 0.1 percent 
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convergence criteria for the iterative analysis and 
found little difference in the optimization results. On 
the other hand, 0.05 percent was the largest con- 
vergence criterion that produced acceptable results 
when global derivatives were estimated by the con- 
ventional approach (finite differencing over the cou- 
pled analysis). A small convergence criterion means 
many additional iterations of the coupled analysis 
and additional computational expense. 

For CSID problems, the GSE approach is effi- 
cient. The coupled analysis is performed only once 
per cycle; thus, the number of structural eigenvalue 
problems to be solved is kept to a minimum. The 
structural code is computationally expensive and 
requires large amounts of computer memory. For 
instance, in the demonstration problem, the finite- 
element analysis required 20 times more central pro- 
cessing time and 20 times more disk space than the 
control analysis required. The GSE approach re- 
quires the calculation of many additional local deriva- 
tives, but most of these calculations involve reduced- 
order control analysis. 

Another strong argument for the GSE approach 
involves human judgment. As a rule, researchers lack 
confidence in global derivatives of a multidisciplinary 
analysis whether the derivatives are calculated by the 
conventional approach or by the GSE approach. On 
the other hand, individual entries in the GSE ma- 
trix (which result from local analyses) generally have 
physical meaning. For example, numerical values 
that indicate the dependence of closed-loop eigen- 
values on open-loop frequencies, or the dependence 
of the inertia properties of the structure on actuator 
mass, may have intuitive appeal. These physically 
meaningful quantities are often available as an out- 
put option from the standard analysis codes. For 
instance, semianalytic formulas for structural deriva- 
tives are common. (See appendix B.) When this is 
not true, the local derivatives can be estimated by 
finite-difference approximation. Finite- difference ap- 
proximations to the control derivatives proved to be 
accurate. (See appendix C.) 

A final argument for GSE, which is mentioned 
in reference 10, is realized in the current research. 
The GSE approach facilitates the smooth operation 
of a multidisciplinary design team. In the current 
research, some team members had more experience 
with control theory, and others had more experience 
with structural analysis. Members used their own 
expertise to develop the local analysis and local 
derivative capabilities. Once individual codes were 
operational, they were combined with generic global 
sensitivity and optimization routines. The end 


result was a powerful and flexible CS1 optimization 
capability. 

Complexity is one definite disadvantage of the 
GSE approach. A large number of local derivatives 
must be communicated to the GSE computer code in 
exactly the correct form and order. Inconsistencies 
within the coupled analysis can be a problem. For 
instance, the structural outputs contain derivatives 
of u; 2 , but the control analysis calculates derivatives 
with respect to a?. 

Another potential disadvantage of the GSE ap- 
proach is the size of the set of linear equations to be 
solved. For the current example, the size never ex- 
ceeded 228 x 228 and thus was not a concern. The 
number 228 represents 2 masses, 6 unique compo- 
nents of the symmetric inertia matrix, 180 elements 
of the 3 mode-slope matrices, 20 elastic modal fre- 
quencies, and 20 unique real parts of the closed-loop 
eigenvalues. The actual size is slightly smaller be- 
cause repeated eigenvalues (due to structural symme- 
try) are eliminated. The size of the GSE matrix can 
expand if the number of modes used in the reduced- 
order model increased or if additional constraints 
(e.g., structural member buckling constraints) are 
added to the optimization problem. 

Comparison of GSE and Conventional 

Approaches 

In the final analysis, the value of the GSE ap- 
proach needs to be assessed relative to the conven- 
tional approach. Clearly, it is easier to develop a 
procedure that perturbs the coupled analysis many 
times and estimates the global derivatives through 
the use of a forward-difference formula than it is to 
calculate and manipulate the local derivatives and 
solve equation (1). However, choosing a perturba- 
tion step size is difficult. The accepted method for 
choosing step size is to test successively smaller steps 
and stop when no significant change in output occurs. 
This method assumes that the coupled analysis is 
perfectly converged; thus, the method requires many 
repetitions of controls and structures codes to satisfy 
a small convergence tolerance. As a consequence, this 
method of determining step size is expensive. More- 
over, the step size that proves to be appropriate for 
a given design point is not guaranteed to remain ap- 
propriate t hroughout the problem domain. 

To compare the GSE and conventional approaches, 
demonstration case 1 was solved by each method. 
As shown in figure 9, the two approaches have sim- 
ilar convergence histories. The computational cost 
per cycle for 15 design variables and 20 modes is 
about the same. The perturbation step size for the 
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conventional approach was set at 0.01 percent of the 
nominal value of each design variable. 

Choosing the proper step size for the conventional 
approach was difficult. Perturbation step sizes that 
were several orders of magnitude larger and smaller 
than 0.01 percent were investigated at the initial de- 
sign point. Each of the three structural design vari- 
ables was perturbed separately. The global deriva- 
tives varied enormously with the perturbation step 
size (i.e., changes in sign and in order of magnitude 
were observed), and the choice of appropriate step 
size was not clear. Eventually, the perturbation step 
size was chosen by comparison of the finite-difference 
derivatives with the global derivatives calculated by 
the GSE approach. This experience suggests that 
the GSE approach is the best choice when the con- 
trol and structural analyses are tightly coupled, that 
is, when an iteration between the analyses is required 
for solution. 

In summary, the GSE approach should be favored 
for problems with a large number of design variables 
and for problems in which the coupling between 
disciplines is strong. When a large number of design 
variables occur, the GSE approach is more efficient, 
and when coupling between disciplines is strong, the 
GSE approach is not sensitive to poor choices of 
convergence tolerance or perturbation stop size. 

Concluding Remarks 

This paper describes the development and imple- 
mentation of a general optimization-based method 
for the design of large space platforms through inte- 
gration of the disciplines of structural dynamics and 
controls. The method is especially appropriate for 
preliminary design problems in which the structural 
and control analyses are tightly coupled. The method 


is significant because it coordinates general-purpose 
structural analysis, multivariable control, and opti- 
mization codes and thus can be adapted to a vari- 
ety of controls-struotures integrated design (CSID) 
projects. The method uses the global sensitivity 
equations (GSE) approach. Although the GSE ap- 
proach has been applied to other preliminary design 
problems, this is the first application to the design of 
a large space structure. 

To demonstrate its capabilities the method is 
used to minimize the total weight of a space plat form 
while maintaining a specified vibration decay rate 
after slewing maneuvers. Although the structural 
model has many simplifying assumptions and the 
number and location of actuators an' fixed, this 
proves to be a challenging design problem. With the 
CSID procedure, the platform is redesigned so that 
the mass distribution and dynamic characteristics of 
the structure enhance the use of rate and position 
feedback by the control system. The CSID method 
must trade stiffness that adds structural weight, for 
control effort that adds weight to the actuators. 
The procedure not only makes a favorable trade of 
structural mass for control effort but also satisfies the 
vibration decay rate constraints. 

This research demonstrates that an integrated 
eontrols-structures optimization method can lead 
to significant mass savings, which may not be re- 
vealed by traditional (single-discipline) design meth- 
ods. The solution of the demonstration problem is 
an important step toward comprehensive preliminary 
design capability for controlled space structures. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
September 27, 1991 
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Appendix A 

Global Sensitivity Equations 

This appendix describes the use of global sensitivity equations (GSE) to calculate derivatives of integrated 
controls-structures response quantities with respect to design variables and emphasizes the application of GSE 
to the demonstration problem. References 10 to 12 discuss the theoretical foundations and other applications 
of GSE. 

In theory, the GSE approach can be applied to coupled problems with any number of contributing analyses. 
The global derivatives are determined by solution of a set of linear equations of the form 


Ah = B (Al) 

where A is a matrix of coefficients, h is the solution vector, and B is a matrix of one or more right-hand-side 
column vectors. 


The size and topology of the GSE are determined by the application. For the demonstration problem 
(eqs. (2)), there are two contributing analyses: controls and structures. The inputs and outputs of these 
analyses are shown in figure 4. There are 15 design variables r and g and thus 15 column vectors in the 
matrix B. There are 42 response quantities M s , M iic t, and A and 206 outputs of the structural analysis that 
influence the control analysis (u>, $ (1) , $ (2) , and J). The preceding counts assume that n = 20. Thus, 

matrix A in equation (Al) has at most 248 x 248 entries. 

The actual number of coefficients that need to be computed is smaller than the maximum number of entries 
in matrix A. The closed-loop eigenvalues A occur in conjugate pairs, and only the real part of each complex 
number is constrained. Generally, no more than 20 of the 40 closed-loop eigenvalues have; unique values or 
unique partial derivatives. Likewise, the open-loop eigenvalues uM are not all unique. Because of the symmetry 
of the structure, some eigenvalues are paired. Furthermore, since none of the design variables or response 
quantities alter the symmetry of the structure, the partial derivatives of eigenvalues with respect to changes 
in the input quantities also occur in pairs. Removing from matrix A row's and columns that correspond to 
nonunique quantities gives a matrix with about 220 x 220 elements. The exact number can change depending 
upon which structural modes are retained in the reduced model. Many of these elements have the value of zero 
and do not need to be computed. 


The topology of matrix A for the reference problem is 

A = 


1(22x22) 


■^ 2 (22xl!)8) 


A 1(I!)8x|) °(198x21) 1(198x198) J 


(A2) 


where the submatrices have the indicated sizes and I is the identity matrix, 0 indicates a submatrix filled with 
zeros, and submatrix A] is defined as 


A'f = 


mtz: 




00 


( 1 ) 




<l h 

A\/a„ • . 


Mi 


(A3) 


where the local derivatives are calculated by a semianalytic method discussed in appendix B. The submatrix 
A 2 is defined as 
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where all local derivatives arc 

approximated 

by finite differencing. 
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Matrix B in equation (Al) has rows that are associated with each unique response quantity and each unique 
output of the structural analysis and are ordered to be consistent with matrix A; the columns in matrix B 
correspond to the design variables. Thus, matrix B has the form 


B = 


where 


Bi 


and 


B, 


°(21x») B 2 ( 2 I X 12 ) 

_ fdlKlxS) °(199xl2) 

r OAL 

0M S 
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- 7)r\ 

<)r> 

ff J 
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oa/ mI 1 
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dAi 

ff 


OX j 
^712 

o\„ 

■ dltn 


On 12 J 

of B, 

art' t 

evaluated 


(AS) 


(AG) 


(A7) 


The solution vector h in equation (Al) contains global derivative's of response and structural output 
quantities with respect to each of the 15 design variables. Thus, for the first design variable (which corresponds 
to the first column vector in matrix B) 


h r = 


dA/a,., 


dr i 


dX\ 

dr\ 


dXn 
'' dr\ 


dM s du)\ 
dr i 


duj n dcj> 


dr i 


1) 

11 


dr i dr] 


dr i 


<Lh 

dr i 


dJ* 

dn 


(AS) 


The- first 2 + n entries in each h, vector contain global derivatives that are needed by the optimization phase 
of the coutrols-structurcs integrated design (CSID) process. 

The solution of equation (Al) is calculated with a standard lower and upper triangular decomposition 

a gorithm. The computer code that solves equation (Al) also estimates the condition number N of the matrix 
denned as 

N = || A (| || A~ 1 1| (A9) 
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where the matrix norm is defined as 


I A|| = max ( Y l 


(A10) 


Condition numbers between 10 :5 and l() fi were typical for the demonstration problems. These exponents suggest 
a loss of 3 to (> significant figures (out of about 14 significant figures on a 64-bit machine) due to numerical 
errors during the solution of equation (Al). The effect of conditioning errors would have been greater than 
this if the individual entries in the GSE matrix were not scaled as suggested in reference 20 tor example, 
dM- u -i/d'h is multiplied by before being placed in the GSE matrix. After equation (A 1) is solved, e 

true value's of the global derivatives can be recovered through reversal of the scaling procedure. 
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Appendix B 
Structural Analysis 

This appendix describes tin* structural model 
used in the demonstration problems. The defini- 
tions of the structural design variables are given, 
and the procedure's used to calculate; the; outputs of 
the structural analysis are explained. The finite- 
element model, (comprised of beam elements and 
point masses, is assembled with the; program Engi- 
neering Analysis Language' (EAL; see ref. 14). 

The geometry of the reference; configuration is 
shown in figure 2. The size and shape of the plat- 
foim does not eliange* during the; optimization pro- 
cess. Structural e;leme;nts arc; divided into three; 
groups (bus. antenna, and antenna support.) and an' 
size'd by the optimization process. The' thm; d('- 
sign variable's /■/ control the' oute'r radius of groups 
ol stiue tuial e'le'ine'nts and assume; value's b('twe;e;n 
i and 10 cm. 


I he truss-like' bus structure is composed of 10 bays 
that are; .4.0 m long by 1 .5 m wide' and high. Figure' 10 
indicates the conne'ction of me'inbe'rs to form a typical 
bay. The; truss ele'inents in the bus arc* tube's with a 
constant, wall thickne;ss of 0.159 cm and with oute'r 
radius controlled by de'sign variable' r\. The bays 
on each end of the bus include; extra me'inbe'rs to 
support the vibration suppression actuators locate'd 
at the centers of the'se two bays. The' actuators are; 
mode'led as point masses t hat, are connected to e'ach 
corner of the' bay. 

The' flat, circular ante'imas with diame'tc'rs of 15.0 
and 7.5 m are; forme'd with 12 radial and 12 circmn- 
fe'n'iitial be'am cle'me'nts. The antenna truss (dements 
have' the same' wall thickne'ss as the bus truss (de- 
ments. but their outer radius is controller! by the de- 
sign v ai indie* 1 2 . The' ante'imas are; supported by truss 
members that connect, the four joints e>n t he* top face 
of the last bay to the center of an antenna. The dis- 
tance' between the centerline of the* bus structure and 
the* plane of e;ach ante;nna is 5 m. The* e>ute;r radius 
of the antenna support memhe'rs is controlled by rp 

4 he* mass and inertia prope;rtie\s of the; re;fe;renc*e; 
configuration and its static response; charactendstics 
are d('te;rmine;d with EAL. Standarel routine's return 
the total mass, the; truss structure mass, and the 
center-of-mass location. The'se* routine's also calcu- 
late* t he inertias about the center of mass. The par- 
tial derivatives of the truss mass and of the elements 
of the; moment, of inertia matrix with re;spcct to r, 
and A/ a( . t are* approximated with a finitevdifferemee; 
technique. 


The dynamic characteristics of the re'fere'iice con- 
figuration are determined through solution of an 
eigenvalue* problem! of the form 


K'P - arM^ = 0 


(Ml) 


where K is the stiffness matrix, M is the mass ma- 
trix, lepiesent t he* eigenvalue's, and 'E repiesemt 
the* eigenveetoi s, normahze'd so that 1VT'E / eejuals 
unity. 


Typical value's of u) 2 arc list cel in t able' 1. The* 
hist 2(i edge*] lvalue's arc* liste'd for structure's that 
correspond to thre;e; different. cycles of de'inoiist rat ion 
e*ase* 1 . The first six eigenvalue's. which are' associated 
with the rigid-body mode's, are* identically zero. Only 
the* eigenvalue's of the* first 20 elastic moeie's are 


transmitted t.o the; controls analysis. Also, the* holel- 
face type' in table* 1 indicates t he* pre'semce* of re*pe'ate*d 
eigenvalue's. 

1 he* paitial derivat ive's of and 'E with re'spe'cf. 
te> the design variable's and t he* mass ed actuators are* 
calculated with semuanalytie formulas de'seribe'd in 


references 18 and 19. For instance', the eigenvalue* 
derivative's with re'spe'cl to rj arc* given by 


Or 1 


V 1 


OK 

— u) 

Or | Or 


2 0 M 
1 J 


1 he* de*i ivat ives e>f the* stillness and mass matri- 


* 

(M2) 

mass 

mat ri- 


ces can be accurately approximated with a finite- 
difference* technique. Although the solution of the 
eigenvalue problem (e*q. (Bl)) is computationally e*x- 
pensive (e*.g., 5 min of CPU time* on a CONVEX 220 
for the* demonstration case*), the partial derivative's of 
eigenvalue's and eigenvectors can be* calculated accu- 
rately and with little e'xpense. cspe'cially if t he* mass 
and stiffness matrices are; updated edliciently. (See 
ie;f. 19 foi implementation derails and cost estimate's 
of various med.hods.) 


Whe'n the solution to expiation (Bl) produce’s re- 
pented e'ige'u value's, then the; calculation of partial 
derivative's of those particular eigenvalue's and of cor- 
responding eigenvectors m pure's speeial care*. (See 
ref. 22.) In the demonstration problem. rc;pe*ated 
eigenvalues are* caused by symmetry in tlie re'fere'iice* 
configuration, riic* mode* shape's associated with re- 
peated eigenvalue's always involve; bending or rotation 
of the antennas. Figure 11 shows the mode shape's 
for a typical pair of repeated eigenvalue's. Since* the* 
bus structure is not involved in either mode, the; ro- 
tational components of the matrix \E at the act uator 
locations must c*qual zero. Moreover, since the de- 
sign variable's do not influence the symmetry of the* 
structure, the partial derivatives of these rotational 
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components with respect to structural parameters 
must also equal zero. Thus, the method developed 
in reference 22 is not required here. 

The outputs of the structural analysis are M s , 
J, u>, and 'I'. The derivatives of these quantities 
with respect to M ac t arc calculated routinely, and 
the derivatives with respect to the structural design 


variables are calculated once per cycle. The matrix d> 
contains modal displacements and rotations at each 
node in the finite-element model. Only the modal 
rotations that correspond to the three actuator loca- 
tions are required by the controls analysis. Therefore, 
these rotations are extracted from the matrix and 
placed in matrices d> 11 ^ . and d» ' ~ 1 . 
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Appendix C 

Actuator Mass Calculation 

The vibration suppression actuators, located at specified structural nodes, are sized so that they damp the 
elastic motion that results from a minimum-time reference rotational maneuver. To determine the actuator 
mass, a reduced-order model of the spacecraft, which is obtained by modal truncation, must first be provided to 
the control analysis from the structural analysis. Next, the vibration suppression control law must be defined. 
A closed-loop model results, which may be put into state-space form. The reference maneuver that excites the 
closed-loop system is next defined, and the response to this excitation is determined. The resulting vibration 
suppression torque histories are scanned to find peak values, and the required actuator mass is determined 
from these values. The remainder of this appendix describes the details of this procedure. 

Eigenvalue analysis of a finite-element model of the spacecraft as described in appendix B allows construction 
of a truncated model in the control analysis. This reduced-order model consists of n elastic modes and has the 
general form 

q + Dq + Aq = T t F + T r T (Cl) 

where q represents an n-vector of unknown modal coordinates, q and q are time derivatives, D is a 
diagonal damping matrix with elements 2 pw,;, p is the modal damping ( p = 0.005), A is a diagonal 
matrix of eigenvalues u >f, and and T r are translational and rotational components of the structural 
eigenvectors at the points where external forces F and moments T are applied. Reference 23 discusses the 
choice of an appropriate number of modes, and reference 24 discusses the classical transformation of the 
differential equations of motion from physical to modal coordinates. The explicit form of the right-hand 
side of this equation is determined once the vibration suppression system and the reference excitation are 
characterized. 

The vibration suppression system employs a robust dissipative control strategy called collocated elastic 
control (CEC), since it provides stability despite unmodeled dynamics and parameter uncertainty. (See refs. 15 
to 17.) Expressions for a single pair of CEC actuators that produce torques T, and T 2 (as in the demonstration 
problem) are developed here, although this procedure may be generalized to any number of actuators. 

Collocated with the CEC torque actuators are attitude and attitude rate sensors. The actual attitude 
at each of the CEC actuator locations includes the rigid-body attitude ao as well as the rotational elastic 
deformation. These angles are 

a l = ao + 3> (1) q and a 2 = ao + <l> ^q (C2) 

Similarly the measured attitude rates are 

aj = ao + <f>^q and a 2 = ao + 4>^q (C3) 

To remove the effect of the rigid-body attitude and rate, the actuators act in a differential mode. Define 


a = ai — a 2 = <3>q 
a = ai - a 2 = 3>q 

where $ = ^l 1 ) - $( 2 ). The actuator outputs are constrained so that Ti = -T 2 = W so that the CEC system 
does not affect rigid-body motion. The control law has the form 

W = — G p a — G r a (C5) 

where the G p and G r are 3x3 positive-definite, symmetric position- and rate-gain matrices. Symmetry of the 
gain matrices is guaranteed implicitly by the following decompositions: 

— y 1 — 'Y 

Gp = G p Gp and G r = G r G r (C6) 
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The 12 unique element s of t he decomposed gain matrices G ?) and G ( are the design variables g. I bus, 


G, 


9 1 

0 

o " 


97 

0 

0 " 

92 

9 :\ 

0 

and Gj- — 

98 

.<79 

0 

_ 9 \ 

95 

%_ 


.910 

.<7ll 

<712 _ 


(C7) 


when' 0 < <1, < 200. 

The closed-loop equations of motion are written with equation (Cl). For now, the reference maneuver is 
given as some moment U applied at a point. The closed-loop elastic equation of motion, which is obtained bv 
considering both the reference maneuver and the CEC system, is 


q + Dq + Aq = # U + <1>^ Ti + T 2 


(C8) 


where the 3 x n mode-slope matrices and <1> (2) contain rotational components of the eigenvectors 

at the point of application of the maneuver moment U and at the two CEC torque actuator locations. From 
$ <p( * 1 — , equation (C8) can be written as 


q + Dq + Aq = U + W 


(CD) 


By substitution for the control law (eq. (C5)) and use of equation (C4), the closed-loop equations of motion 

become 7 

q + Dq + Aq — 4> |0) U 


(CIO) 


where' positive definite matrices D and A are 


D = 
A 


I) + # r G r # 1 
A + * r G p * j 


This system may easily be written in state-space form as 

x = Ax + BU 

for the 2 n states x = (q r ,q 7 ) 7 '. The 2 n x 2 n matrix A is 


0_ I_ 
-A -D 


(Cll) 


(C12) 


(Cl 3) 


B 


(C14) 


The 2 n x 3 matrix B is given bv 

0 

$(0) r 

The closed-loop eigenvalues, which are used as constraints in the optimization process, are the eigenvalues of 
matrix A. 

The response x(*) of the closed-loop system to a reference excitation dictates the required vibration 
suppression torques. This reference excitation was chosen to be a bang- bang type, minimum-time lotational 
maneuver given by 

Ui 1 f u max (0 < t < tf/2) 

U(t) = U2 — S —U m ax (f//2 < t < t f) 

_U 3 J [o ( t>tf ) 

where' Umax is the maximum applied maneuver moment given by 

4 


(C15) 


Umax — m a f 


(016) 
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Here, Uj represents components of the maximum maneuver moment, t f is the maneuver time, J is the spacecraft 
inertia matrix, and acj is the vector of (small) desired attitude angles. (The reference excitation, which depends 
on J, is different for each repetition of the controls analysis.) The initial orientation angles and times are 
assumed to be zero. Maneuver parameters tj and ctj remain fixed throughout the optimization procedure. 
The reference maneuver only exists to simulate a realistic disturbance to the spacecraft. Thus, the mass of the 
maneuver actuator (assumed to be zero) is never changed. 

The actuator mass is determined by solving equation (C12) for the response x(t) and then calculating from 
equations (C4) and (C5) the actuator torque W(f), which is a function of x. The response is given by 


x(0 = cxp(At)xo + ex 


P(A0 I' 
J o 


exp( — At)BU(t) dr 


(Cl 7) 


Evaluating the convolution integral yields 


A 1 [ exp(Af) - I] Bu tnax (0 < t < tj/2) } 

x(t) = <j A 1 {ejcp(At) - 2 exp [A(/. - t f /2)\ +1} Bu Iliax (t f / 2 < t < f f ) 


A' 1 exp(At) [exp(— Aty/2) - I] 2 Bu max 




(C 18 ) 


This expression may be used to find the torque at any time by recalling that the C’EC control law (eq. (C5)) 
may be rewritten as 

W(0 - | - G r ^] x{t) (Cl 9) 

For the demonstration problem, the peak values of W are found simply by scanning each required control 
torque history for each of the three coordinate directions for a sufficient period of time after the maneuver 
ends. Limiting the search to times greater than tf allows calculation of the unforced response recursively at 
minimum computational cost, (c.g., 10 sec of CPU time on the CONVEX 220). Thus. 


X/+ 1 = «xp(A At.)x; (/ - 0, 1, 2, n t ) (C20) 

where x 0 = x{tj), and At is no greater than P/2, where P is the period of the highest vibration frequency in 
the truncated model, and u/ is the number of points scanned. For the demonstration eases, At — P/\iS and 


the time histories are scanned starting at tj for a time period equal to twice t he period of the lowest, vibration 
frequency. 

Having identified the peak torques, the actuator mass is determined with a functional or empirical 
relationship between torque output and actuator mass. For the demonstration problem, a constant of 
proportionality is assumed, and the contributions to the total actuator mass from each coordinate direction 
are 

rn-if = 2M r {w\,w$, w$) T = 2AI r W* (C21) 

where ie ; pwik torque occurring at time t* for the 7t.l1 coordinate direction, and A It is 

tlu 1 constant of proportionality hetween mass and torque output (A/7 1 = 1.0 lbm/ft-lbf = 0.304 kg/N-111). The 
factor of 2 exists because there are two CEC actuators. The actuator mass then is 


A/act = mi + tri i 2 + m 3 


(C22) 
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Appendix D 

Computer Implementation 

The integrated controls-structurcs optimization 
procedure is implemented as a batch job on a 
CONVEX 220 computer. The procedure consists of 
a set of FORTRAN codes that are linked by oper- 
ating system commands and by input and output 
(I/O) files. The operating system on the CONVEX 
computer is a derivative of UNIX, and the command 
language interpreter is similar to the C shell (csh) 
available on many machines. Table 2 contains a list 
of all FORTR AN codes and their I/O files. 

The batch job developed for the current research 
is useful as a model of a generic optimization pro- 
cedure for coupled, multidisciplinary optimization 
problems. This batch job can be modified to ad- 
dress a variety of controls-structurcs integrated de- 
sign (CSID) optimization problems. Figure 12 is a 
list of the batch submittal file, and figure 13 is a 
conceptual flowchart that explains the operation of 
the batch job. In figure 13, the rectangles represent 
FORTRAN programs, the diamonds represent logical 
tests, the heavy lines indicate normal data flow, and 
the dashed lines indicate a shortcut, which is used 
whenever possible. The numbers in parentheses be- 
side the flowchart correspond to numbered comments 
in figure 12. 

The command language implementation scheme 
has a number of advantages over the traditional 
programming language (e.g., FORTRAN) implemen- 
tation. This is a time-consuming procedure (e.g., 


one demonstration problem required about 25 hr of 
CPU time on the CONVEX computer). This com- 
mand language implementation facilitates monitor- 
ing, restarting, and debugging the process. Since 
each execution of each FORTRAN program creates 
new output files, a simple list of file names and times 
of creation is used to deduce the state of an execut- 
ing process. At the same time, any output file can be 
examined or printed to obtain information about the 
progress of the optimization. If the optimizer or anal- 
ysis routines do not appear to be operating correctly, 
the process can be stopped and restarted at the cur- 
rent step or at a previous step. To restart at the 
beginning of any cycle, only the files “VARIABLE,” 
“M.act,” and “history.old” need to be edited to con- 
tain the design point, the initial estimate of actuator 
mass, and the history of previous optimization cy- 
cles. If the process stops abnormally, then files such 
as the Engineers Analysis Language (EAL) input and 
output files, which normally are not printed, are ex- 
amined for the cause of the failure. 

A final advantage of this command language im- 
plementation is computer resource management. As 
each FORTRAN execution is completed, the batch 
job condenses the output and deletes as many files as 
possible. Thus, complete information is available for 
debugging when the process stops abnormally, but 
the disk space is saved when the process executes 
normally. The filtering tools available in a UNIX- 
based operating environment (e.g., “grep”) are used 
to extract a summary of information in each output 
file before the file is deleted. (See fig. 12.) 
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Table 1. Eigenvalues Calculated During Optimization Cycle 
in Demonstration Case 1 


These cycles represent the initial and final conditions 
and one intermediate result; bold-face type indicates 
repeated eigenvalues 



Eigenvalues for — 

Mode no. 

Cycle 0 

Cycle 6 

Cycle 27 

1 

0 

0 

0 

2 

0 

0 

0 

3 

0 

0 

0 

4 

0 

0 

0 

5 

0 

0 

0 

6 

0 

0 

0 

7 

17.6 

34.0 

18.2 

8 

37.7 

123.1 

36.6 

9 

71.5 

195.9 

75.5 

10 

86.7 

272.9 

85.1 

11 

86.7 

310.4 

85.1 

12 

89.0 

310.4 

101.8 

13 

156.8 

333.6 

159.4 

14 

407.5 

661.3 

325.6 

15 

407.5 

1524.0 

399.1 

16 

442.0 

1531.0 

399.1 

17 

530.1 

1531.0 

506.1 

18 

787.6 

1914.0 

821.0 

19 

1114.0 

2797.0 

1114.0 

20 

1350.0 

4954.0 

1320.0 

21 

1350.0 

4954.0 

1320.0 

22 

1387.0 

5216.0 

1360.0 

23 

1387.0 

5264.0 

1360.0 

24 

2588.0 

5264.0 

1871.0 

25 

2588.0 

6733.0 

2061.0 

26 

2600.0 

9102.0 

2531.0 
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Tabic 2. FORTR AN Programs With Corresponding Input and Output 
Files Used by CSID Procedure 


Program name 

Input files 

Output files 

preeal 

M.act, modes. old, 
VARIABLE, eal.old 

modes.dat, eal.in 

EAL version 324 

cal. in 

eal.out 

posteal 

cal. out 

modes.dat, gse.dat 

MACT 

modes . da t , act co n . < la t , 
VARIABLE 

actmass.out. M.act 

dMACT 

modes . d at , ae t eon . d at 
history.old. VARIABLE 

history. new, cderv.dat 

gsebld 

gse.dat, cderv.dat, 
actcon . dat , VA R I A B L E 
rnodes.dat 

DESVARBL 

optimz 

DESVARBL, FIXEDPAR 
history, new 

VARIABLE, M.act 





Cycle 



(a) Conventional CSID procedure, (b) Alternate CSID procedure that uses global 

sensitivity equations. 

Figure 1. Optimization approaches for coupled controls-structures integrated design (CSID) problems. 


• 3-axis torque actuator 



Figure 2. Reference configuration of geostationary platform with three-axis torque actuators. 
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Figure 3. Flowchart of iterative procedure for coupled analysis. 


r, M act 

©, J. g 

* 


Structures 

Controls 

* 


M s , co, O, J 

X, M ac( 


Figure 4. Input and output quantities for contributing analyses. 
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(b) Optimization history of inaximnm closed-loop eigenvalue compared to its allowable value. Constraint 
is initially violated. 


Figure 5. Results of demonstration case' 1. 












(b) Real part of maximum closed-loop eigenvalue. 

Figure 8. Comparison of actual computed values of selected outputs of coupled controls-structures analysis 
with linear approximation to those values. 
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Figure 9. Comparison of optimization histories of demonstration ease 1 with GSE and conventional 
approaches. 



Figure 10. Typical hay of truss structure for reference configuration. 



• Actuator locations 




Figure 11. Typical mode shapes for eigenvectors of repeated eigenvalues. 



#! /bin/csh 
#6$-o output -eo -nr 
# 

# INITIALIZATION 

# 

# copy input files to scratch disk 

# 

date 

if ( -f actcon.dat ) then 

echo RESTART do not copy anything 
else 

echo copy all files needed for optimization 
cp -/eal.new eal.old 
cp ~/M.act M.act 
cp -/VARIABLE VARIABLE 
cp -/FIXEDPAR FIXEDPAR 
cp -/history . old history. old 
cp -/preeal preeal 
cp -/posteal posteal 
cp ~/MACT MACT 
cp - / dMACT dMACT 
cp -/gsebld gsebld 
cp -/optimz optimz 
cp -/actcon.dat actcon.dat 
endif 
rm L* 

rm modes . dat 
@ numopt = 0 
# 

# (9) new cycle? 

# Repeat the outer loop as long as a file "VARIABLE" exists 

# 

while ( -e VARIABLE ) 

# 

# reset counters 

# 

rm DESVARBL 
@ cycle = 0 
@ ealcnt = 0 
0 macnt = 0 

# 

# (5) converged? 

# 

# Iterate between EAL and MACT until new mass of actuator 

# and previous mass of actuator are close enough 

# Repeat inner loop as long as file "M.act" exists 

Figure 12. Batch submittal file for implementation of CSID procedure with UNIX command language. 
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# note: inner loop terminates abnormally after 7 calls to EAL 

# or after 50 approximations to EAL output 

# 

while ( -e M . act ) 

@ cycle++ 

echo CYCLE $cycle 
if ( $ealcnt >= 7 ) then 
echo MACT output 
cat actmass.out 
exit 1 

endif 

if ( $macnt >= 50 ) then 
echo actmass output 
cat actmass.out 
exit 1 

endif 

# 

# (1) EAL PREPROCESSOR 

# 

preeal 

echo preeal completed 
mv M . act mass. act 

# 

if ( -e modes.dat ) then 

# 

# ACTUATOR MASS CALCULATION 

# use estimated frequencies 

# 

@ macnt++ 

echo start MACT approximate analysis $macnt 
MACT 

echo MACT completed 
grep "mass" actmass.out 
mv modes.dat modes. old 

else 

# 

# (2) EAL Structural Analysis 

# calculate actual frequencies 

# 

@ ealcnt++ 

echo start EAL full analysis $ealcnt 
~eal/eal324/eal324 <eal . in >eal.out 
echo EAL completed 
rm L* 


Figure 12. Continued. 
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# 

# (3) POSTPROCESSOR 

# 

cat eal.outlgrep -v DATA I grep -v EX >eal.log 

echo Prepare input to actmass 

posteal 

echo posteal completed 

# 

# (4) ACTUATOR MASS CALCULATION 

# 

echo Calculate actuator mass 
MACT 

echo MACT completed 
grep "mass" actmass. out 
mv modes.dat modes. old 

endif 

end 

echo converged 

echo DESIGN VARIABLES 

cat VARIABLE 

echo Structures model 

grep " J=" eal.in 

echo Structures output 

cat modes. old 

# 

# (5) END INNER LOOP 

# 

echo 

# 

# (6) dMACT needs modes.dat to estimate f.d. derivatives 

# 

cp modes. old modes.dat 
dMACT 

echo Controls output 
cat actmass. out 
echo dMACT completed 

# 

# (7) gsebld assembles and solves GSE 

# 

echo output from gsebld 
gsebld 

echo gsebld completed 


Figure 12. Continued. 
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=tt= =tt= =#= =#= 4^ =#= =4*= =#= 


(8) Global derivatives are in a file named "DESVARBL" 

This file not produced if design is feasible and objective 
stops improving 

@ numopt++ 

echo BEGIN OPTIMIZATION $numopt 
if ( -e DESVARBL ) then 
cat DESVARBL 
rm VARIABLE 

else 

exit 1 

endif 

(9) check maximum number of optimization cycles 

if ( $numopt >= 10 ) then 
opt imz 

echo FINAL CYCLE optimz completed 
echo output from optimz 
cat OPTRESLT 

echo HISTORY of OPTIMIZATION 
cat history. new 
mv history. new history. old 
echo completed 
exit 1 
else 

optimz 

mv history. new history. old 
echo optimz completed 
echo output from optimz 
cat OPTRESLT 

# 

endif 

# 

# if modes.dat exists - the inner loop will skip EAL 
rm modes . dat 

# 

date 

end 

# END OUTER LOOP 

# 


Figure 12. Concluded. 
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Figure 13. Flowchart of CSID procedure with GSE approach. See table 2 for FORTRAN programs and 
figure 12 for numbered comments. 
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